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(54) Adaptive noise cancellation device 

(57) An adaptive noise cancellation device comprises: convolution logic 10 for convolving the signal from a 
reference input x with a discretised L-tap filter to form a filtered reference signal; and logic 20 for subtracting 
the filtered reference signal from a signal input to form an output signal; logic for generating the filter taps as 
a linear combination of N basis functions each having a corresponding coefficient Ck; and logic for repeatedly 
determining the coefficients Ck which minimise the power in the output signal z, characterised in that N is less 
than the number of filter taps L and the basis functions have a portion of finite width, outside of which portion 
the functions are substantially zero, both in the frequency and time domains, in an embodiment they are 
gaussians. A full-duplex speakerphone is disclosed including such a noise cancellation device. 
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ADAPTIVE MOISE CANCELLATION DEVICE 



5 The invention relates to adaptive noise cancellation devices and to 

the use of adaptive noise concellation devices as howling cancellers for 
full duplex speakerphones . 

It is common practice to use adaptive noise cancellation devices in 
10 the transmission of intelligible speech from an -environment with 

considerable acoustic noise, such as an aircraft cockpit, a car etc., 
utilising as a reference signal the signal from an additional microphone, 
located sufficiently far from the speaker to be insensitive to the speecli 
of the speaker and as close as possible to the noise source. Adaptive 
15 noise cancellation devices are also used in full-duplex speakerphone 

applications to prevent howling resulting from the feedback path between 
speaker and microphone of the speakerphone device. 

Various techniques for adaptive noise cancelling are known, some of 
20 these are described in a paper by B Widrow et al, "Adaptive Noise 

Cancelling - Principles and Application of the LMS adaptive filter", Proc 
IEEE, Vol 63, No 12, Dec 1975, pp 1692-1716, and in a paper by II R 
Sambur, "Adaptive Noise Cancelling for Speech Signals", IEEE trans. ASSP, 
vol 26, 1977, pp 419-423. 

25 

An adaptive noise canceller first produces an estimate of the 
characteristics of the transf ormal:ion from the reference noise signal to 
the noise component of the signal at the main microphone. This 
transformation actually depends on two acoustic coupling paths: the 
30 first between the noise source to the additional microphone, the second 
between the noise source to the main microphone. Then the reference . 
noise signal is used to model the noisy signal component at the main 
microphone. This noise component is subtracted from the actual 
microphone signal. In the absence of the speaker signal and assuming 
35 that the tranrf ormatio:* is exactly identified, the difference between the 
actual main microphone input and the estimated noise at the adaptive 
filter output would be zero. In the presence of the speaker signal this 
difference signal contains mainly the signal from the speaker - 

40 The filter characteristics are dynamically changed for optimum 
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removal of the noise signal. The two acoustic paths, which define the 
transformation between the additional and the main microphone are 
generally nonstationary - The adaptation rate has to be sufficiently high 
to adaptively track changes in the impulse responses of these paths due 
5 to the movement of people or objects in the room, the movement of the 
noise source or the change in the noise characteristics - 

Adaptive filters usually employ the well established "LMS 
algorithm", which is known also by the name "stochastic gradients", see, 
10 for example, Widrow B, et al "Stationary and Nonstationary Learning 

Characteristics of the LMS Adaptive Filter", Proc IEEE, August 1976, Vol 
64, No 8, PP 1151-1162. 

According to the LMS algorithm, the output of the adaptive filter 
15 is required to be as small as possible in the sense of least-mean-square 
error, ie the output power is minimised. The LMS algorithm updates N 
unknown filter coefficients each input sample and produces an 
approximately optimal solution in 0(N) computations. The LMS algorithm 
works best when the reference signal is a white noise signal. However, 
20 for actual signals which differ from a white noise, and are 

non-stationcjry , the convergence of the LMS method is very poor. 

The required filter length, ie the number of taps N, is determined 
by the length of impulse responses of the acoustic leakage paths, that is 

25 by the reverberation times of the room. Typical reverberation time of 
most rooms is less than 400 msec. The adaptive filter length must 
therefore be ;>bout 100 to 200 msec, ie in the range of 1000 to 2000 taps 
for a common rhoice of sampling frequency of 8 KHz. Not counting the 
filtering operation itself, the computational requirements for adaptive 

30 LMS filter implementation are therefore at least in the range of 16 to 32 
million instructions per second - a very substantial amount of computing 
power. 

Another approach to the problem of adaptive noise cancellation is 
35 known as the sub-band acoustic echo canceller (SBAEC) , see for example a 
paper by Andre Gilloire, "Experiments with Sub-Band Acoustic Echo 
Cancellers for Teleconferencing", Proc of IEEE ICASSP 87, April 1987, pp 
2141-2144, These have advantages over the standard LMS method both in 
adaptation rate and in computational complexity. They do, however, have 
40 a residual error due to channel interdependence, which has not been taken 
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into account in the adaptation scheme, and this limits the number of 
channels and therefore the computational improvement which may be 
obtained. 

5 This invention provides an adaptive noise cancellation device having 

a signal input, a reference input and a signal output, the device 
comprising: convolution logic for convolving the signal from the 
reference input with a discretised L-tap filter to form a filtered 
reference signal; and logic for subtracting the filtered reference signal 
10 ' from the signal input to form the output signal; logic for generating the 
filter taps as a linear combination of N basis functions each having a 
corresponding coefficient Ck; and logic for repeatedly determining the 
coefficients Ck which minimise the power in the output signal, 
characterised in that N is less than the number of filter taps L and t!ie 
15 basis functions have a portion of finite width, outside of which portion 
the functions are substantially zero, both in the frequency and time 
domains . 

The LMS formulation forms the basis for the approach of this 
invention- Iio;7ever, unlike the prior art, the filter is represented as a 
sum of N basis functions, where N is less than the filter length. This 
is a generalisation of the SBAEC, In the latter, the signal is 
prefiltered by several band filters. However, for each sub-band a 
conventional LMS algorithm is applied to adapt a number of filter 
coefficients in each band, the interdependence of the different bands 
being neglected. By contrast, in the approach of this invention there 
are as many "bands" as degrees of freedom. In each sub-band there is 
only one filter coefficient to be defined and interdependence of adjacent 
bands is syst^jmatically accounted for. 

Preferably the basis functions have a small time-bandwidth product. 
That means that the basis functions have to be well concentrated in the 
time domain, ie their time-width should be as small as possible and 
simultaneously these functions have to have a bandwidth in the frequency 
domain as small as possible. 

In practice what this means is that the basis functions are such 
that within the desired accuracy, both the basis functions and their 
Fourier transforms may be approximated by narrow functions with finite 
40 supports. The use of such basis functions will mean that the filter 



BNSOOCID: <G8 2274372A_L> 



25 



30 



which is obtained by combining them, will in general be non-causal, 
resulting in some additional delay. However, by a judicious choice of 
parameters this delay can be made equal to only a small portion of the 
filter length, and small compared to other delays in the system. 

5 

In one embodiment, the basis functions are generated by an equation 
of the following form, 

u'jt(t)=v(t)exp{jA-(Aw) t) (1) 

10 such as are used in a discrete Fourier transform. Here, j = ^T-l, k is 
an index which runs from -N/2 to (N/2)-l, and A<a is the step in the 
frequency domain. In this form the requirement of small time-bandv;idth 
product is applied actually only to the window function w(t). 

15 The requirement that the basis functions have small time bandwidth 

product results in negligible interdependence between basis functions 
with indices i and k where |i-k| is larger than about 3, Therefore the 
system of equations for filter coefficients, which is obtained by a 
differentiation of the expression for mean-square error, has a sparse 

20 matrix, in which only diagonals near the main diagonal differ from zero. 
Such system can be solved efficiently with 0(N) computations- Since a 
computation of the system matrix itself requires a calculation of 
correlation functions, which is performed in O(NlogN) operations, the 
whole filter adaptation can be made in O(NlogN) operations per 

25 computational block, yielding a new accurate update of the filter 

coefficients per block, assuming that the block length is some small 
integer multiple of the filter length. 

An embodiment of the invention will now be described, by way of 
30 example only, with reference to the accompanying drawings, wherein 

Figv 1- is a schematic block diagram of the adaptive noise 
cancellation device of the embodiment of the invention; 

35 Fig. 2. is a schematic block diagram of a speakerphone which 

includes an adaptive noise cancellation device as a howling canceller. 
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A block diagram of the adaptive noise cancellation device is sliown 
in Figure 1. The device accepts as input the audio signals reference 
input signal x and main input signal y, and produces as output the signal 
z. The device includes the following elements: the convolution logic 10, 
5 which convolves the reference signal x with the current filter 

coefficients; the subtractor 20, which. subtracts the convolver output 
from the main input signal y (appropriately delayed); the correlation 
logic 30, which computes the autocorrelation function of the reference 
signal x and the cross-correlation function between the reference and the 

10 main input signal; the equation set logic 40, which computes, from these 
correlation functions, the coefficients of the system of linear equations 
for filter coefficients; the equation solver 50, which solves the system 
of linear equations; and the control logic 60, which computes short time 
pov;er or amplitude estimates for the reference and the subtractor output 

15 signals correspondingly and determines a control parameter for the 

correlation block. Finally, the automatic gain control (AGC) block 70 
multiplies the subtractor output by slowly changing gain, the output of 
the AGC block 70 being the output of the device. 

20 In this embodiment of the invention the logical elements are 

implemented by a suitably programmed combination of an IBM Personal 
System/2 personal computer, a plug- in accelerator card which includes an 
i860 processor chip manufactured by Intel, and an IBM Audio Capture and 
Playback card for inputting and outputting the signals (IBM and Personal 

25 System/2 are trademarks of International Business Machines Corporation). 
Using this arrangement it has proved possible to attain real-time 
performance with filters of length 1024. 

However, the invention could equally be implemented by a suitably 
30 programmed general purpose Digital Signal Processor which performs real 
time processing of incoming digital signals, or by using special purpose 
hardware. It will be understood that the blocks indicated schematically 
in Figure 1 and the above description could in practice be intermixed. 
For example, the convolution logic could be implemented in the frequency 
35 domain and rould use the same Fourier transform data and algorithms, as 
those which are used by the correlation logic. 



The reference signal x is fed to the device through a suitable 
analog/ digital (A/D) converter. The input signal y is connected through 
40 another A/D converter to the device. 
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The convolution logic 10 generates a convolution of the reference 
signal, x, with a filter response function Q which simulates the acoustic 
path betv;een the reference noise source and the main signal input. The 
output of the convolution logic 10 simulates the interfering noise signal 
entering the microphone. 

According to the LMS (Least-mean-square) formulation, which forms 
the basis of the present approach, the output of the adaptive filter is 
required to be as small as possible in the sense of the least-mean-square 
error. Since the interfering noise and the desired signal are assumed to 
be uncorrelated signals, the minimisation problem is to determine Q such 
that {ly-QxI}, ie the output pov;er, is minimised, where y is the primary 
signal at the main microphone and Qx is the output of the filter Q given 
the input x to the filter. The filter output Qx = q'^x, where «> denotes a 
convolution operation and q(t) is a filter impulse response. 

The filter Q is represented as a sum of M basis functions, with 
unknown coefficients, Cj^^ The filter has the formi 



where aw is a frequency step, t is time and w(t) is a window function. 

In the frequency domain the Fourier transforms, fj^, of the basis 
functions, <^i^(t), are shifted copies of the Fourier transform W(6>) of the 
v/indow function w(t) ; 



<t jtC t) =v{fc)exp(jic(A€o) t), 



(2) 



(3) 



Hence, the filter frequency response is found to be a combination of 
shifted window responses : 



7 



^;i(o)) = 52 C^,W(u^-kc(o) (4) 



2 



The solution of the minimisation problem leads to a set of M linear 
equations for Cj^: 



J: a,,C,-r, i - -4 (5) 



5 

v;ith the system coefficients aj;^, and the right-hand side ry. If a^]^ and r 
can be determined, then can be calculated. 

In the general case, the system coefficients ajj^ can be express?}d 
10 as follows; 



= J A" (oj) A't w) W (cij-iACij) Pi^(gj-/cag)) <io) , (7) 



15 

where X(6)) and y(o) are the Fourier transforms of x and y respectively 
and the right-hand-side coefficients rj can be expressed as follows. 



r. r-. I <&V(«)^'(w) i'{o})dG> (8) 

20 

(9) 
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By denoting a combination of window responses as : 

K(iOr\x) = f^^* (o>-p Aa>) ft^(a> + iiACo) (10) 
and correspondingly an auxiliary function in the time domain: 

;-(f:,^) = (11) 
the coefficients a^i^ and rj can be transformed: 

^ik = ^^xx^^^^^^'-^^) G^(o>,\i) ^y[R^(x)k(r.ii)] (12) 

and 

= G^CiAOj) = Sr[R^^{z)ur{x)] , (13) 



where R;^;^ is the auto-correlation of the reference signal x and R^y is a 
15 function of the cross-correlation between the microphone output y and the 
reference x. 

Thus ajjf and r^ and the unknown coefficients, Ci^, in the filter 
response function can be determined from auto-correlations and cross- 
20 correlations of the input signals. The filter can then be computed by 
solving equation (5) which minimises the noise interference in the 
output, z. 

The computation of the filter coefficients involves the correlation 
25 logic 30, the equation set logic 40 and the equation solver 50. The 
correlation block 30 updates current auto- and cross-correlation 
functions for further analysis, the equation set block 40 uses the 
updated correlations for setting coefficients of the system of equations, 
and the equation solver block 50 computes the unknown coefficients. 

30 

The window function, w(t), which is used to generate the basis 



5 



10 
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functions in equation (2) is required to have a small time-bandwidth 
product. This results in negligible interdependence between basis 
functions with indices i and k, where |i-k| is larger than about 3, v;her 
i and k are complex exponent indices. Thus the system of equations for 
filter coefficients has a sparse matrix, in which only several, D, 
diagonals above and below the main diagonal differ from zero. Such a 
system can be solved with 0(N) computations. Since a computation of the 
system matrix itself requires a calculation of correlation functions, 
which is performed in O(MlogN) operations, the whole filter adaptation 
can be made in O(NlogN) operations per computational block, yielding a 
new accurate update of the filter coefficients per block. Generally the 
block length should be some small multiple of the filter length, L and 
the filter length, L, is of the order of a few N, where N is the number 
of degrees of freedom. 

The correlation block subdivides the two signals: the reference x 
and the main input y, into frames of constant length that are equal to 
the filter length, L. Odd segments are padded by L zeros on the left and 
even segments are padded by L zeros on the right. 

The Fourier transforms F^^, F^y and Fyj^ of the three one-sided 
correlations R^, R^y and Ry^ all with length L, the lag m being in the 
range 0<m<L*l. are initialised to zero. The Fourier transforms cf the 
padded segments are then used to update the Fourier transforms F^x* Fyy 
and Fyx as follows; 

= «^^xy + ^ xryj. (14) 

v/here Xj.j and Yj,] are the Fast Fourier Transforms (FFTs) of the previous 
segment and Xj and Yj are the FFTs of the current segment. The 
"forgetting factor", a, is a real coefficient in the range 0<a<l . Its 
value is determined by the control logic 60. 

The output of the correlation logic are three Fourier transforms 
F^, Fxy and Fy;^, each containing 2L complex numbers. 



The equation set block is activated every few frames. It consists 
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of a lag domain correlator, a matrix computation block and a right-hand- 
side computation block. 

The lag domain correlator computes the auto-correlation R^x and the 
cross-correlation R^y in the lag domain (- (L-1 )<0<L-1 ) by using the IFFT 
(Inverse Fa^t Fourier Transform) on the outputs F^^. F^y and Fy^ of the 
correlation block. The first L numbers out of each 2L are kept. The 
results are arranged into the auto-correlation, R^, whose values at 
negative lags are copied from values of positive lags, since is 
symmetric and the cross-correlation, R^y, whose values at positive lags 
are taken as the first L numbers of IFFT of Fj^y and values at negative 
lags are taken as the first L numbers of IFFT of Fy^. 

The cross-correlation, R^y, in the range -(L-1) to L-1, are fed 
into the right-hand-side computation block. The computation of the rj is 
based on equation 13 as follows. 

The cross-correlation data are multiplied by a sampled window 
function, and are then padded by zeros to obtain 2L numbers. They are 
then ordered as follows: Rxy(O) , Rxy(l ) Rxy(L-l) ,0 ,Rxy(- (L-1) ) , 

The data are transformed to the frequency domain using FFT and the 
results are decimated by (L/N):l, producing N complex numbers r^. 

The auto-correlation, R^, in the range -(L-1) to L-1, are fed into 
the matrix computation block. The computation of the a|j is based on 
equation 12 as follows. 

The auto-correlation data are multiplied by precomputed samples of 
one of several auxiliary functions, k(T,^), which are computed from the 
window function, for ^1=0 , 1/2 , 1 , . . . ,D/2 , where D is the number of 
diagonals above (and correspondingly below) the main diagonal which 
differ from zero. The window function w(t) and the frequency step aw 
could be chosen so that D was small. 

The D results are padded by a zero in order to obtain 2L numbers 
for each value of ^• They are then ordered as in the Rxy case. The data 
are transformed to the frequency domain using FFT for each 



The FFT results are decimated by (L/2M):1, yielding 2N real numbers 
for each |a. The results for |a=0 are used to fill the main diagonal of the 
system matrix. The result for M=l/2 are used to fill the diagonal just 
above the main one. Complex conjugates of these numbers are used to fill 
the diagonal just below the main one. The results for the other values of 
\x are used to fill the next diagonals, until D diagonals above the main 
one and D diagonals below the main one have been filled. 

The resulting sparse system could be solved by using complex 
arithmetic by Gauss elimination in 0{hl.D^) complex operations. In the 
general case, however, matrix normalisation is performed in the following 
steps; 

1. Matrix symmetrization : The complex system matrix [aji^] is 
multiplied from the left on the matrix [a^^]. D additional diagonals will 
all have small entries of the order of magnitude of those entries which 
were neglected in the main system. They are thus neglected. So in the 
resulting real symmetric matrix only 2D+1 diagonals are left. 

2. Right-hand-side recomputation : The right-hand-side is multiplied on 
the matrix [a*i^j]. 

3. Matrix Regularisation : Theoretically the matrix is always positive 
definite • However this property can vanish due to truncation errors in 
computationr. . Thus it is worthwhile to slightly change the main diagonal 
entries in order to provide positive definiteness of the system matrix. 
The maximum number of the main diagonal is multiplied by a small factor, 
e, and the result is added to all numbers on the main diagonal. 

The above description of the general case has been included to 
illustrate how the device would work for a general window function. 
However, it is possible to choose a window function which simplifies the 
computation. An example of such a function is the gaussian window 
function , 

;^(0 = exp{-((t-C„)x^)^) t^=Illh (IS) 



which the inventors use in the preferred embodiment. This is an 



12 

advantageous choice for the window function because it has the minimal 
possible time-bandwidth product (if that product is defined as a product 
of the second moments of w(t) and W(w)). Here fg is a sampling frequency, 
Y is a parameter, which is chosen so that the practical support of w(t) 
5 is finite; for example, if Y=fs( t2"-t2 )/6 , then v/(t)<0.0002 for t<t[ and 
t>t2, while w(to)=l. The centre of the window function is denoted by tg. 

The Fourier transform of this v/indow function is; 
f>/(«) = yWY^^xp{-.^^!iiii}exp(a>tJ (16) 

4 

10 

The Fourier transforms 'fj;(w) of the basis functions 0k(t) are 
shifted gaussians. 

When the number, N, of degrees of freedom for the filter increases, 
the number of gaussians $(«) on the frequency axis and their density 
increase as well. With the correct choice of parameters (y and a^j) the 
overlapping between gaussians is small. This means that must be of the 
order of 1/N and y must be of the order of N. 

The set of parameters used in the embodiment is ao = 2TTfs/N and y 
in the range from 0.25N to 0,75N. The optimal value for y has been found 
to be Y = (15/32)N, Hov/ever other sets of parameters could be chosen, 

25 With such a choice of parameters, only the diagonals which are near 

the main diagonal in the matrix [ajj^] will have significant values which 
need to be taken into account in the solution. Typically there are three 
non-zero diagonals from each side of the main one. The matrix is sparse 
and can be inverted in 0(N) operations, 

30 

The general formula for the matrix coefficients a^i^ yields for a 
gaussian window function; 

3^;,= Cexp{-i^l(A<o)2{i-;c)^>G^.(-^Aw) 

o 2 



15 



20 



(17) 
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and the formula for the right-hcind-side coefficients rj yields for a 
gaussian window function; 

v/here C=y 77/2. 

The Fourier transforms are sampled in steps of aw in (19) or am/2 
in (18). 

In the case of the gaussian window the steps of matrix 
symmetrisation and right-hand-side recomputation are unnecessary, and so 
the calculacion of the system matrix is much simpler. Calculation of the 
matrix [ajj^] is based on equation 17 as follows; 

The auto-correlation data, Pyy/ in the range -(L-l) to I.-l are 
multiplied by the sampled gaussian and padded by a zero to obtain 2L 
numbers. The data are ordered as in the general case. The data are 
transformed to the frequency domain using FFT, Due to the symmetry of the 
auto-correlation the result will be real. 

The FFT result is decimated by two, yielding 2N real numbers: 
ForF] , . . . ,F2H-i , each of which corresponds to a certain value of ink and is 
placed on one of the 2N secondary diagonals of the matrix [aj].]. The 
multiplication by another gaussian, Cexp{-a (i-k) ^ } , a=(YAw)2/8, is 
actually performed along each of the secondary diagonals. For each of the 
principal diagonals, the value of |i-k| and the value of the gaussian 
multiplier are constant. The maximum of the gaussian is located on the 
main diagonal. Due to a practically finite support of the gaussian only a 
few diagonals are considered as non-zero. 

The matrix is defined on a torus, with the right-hand-side 
connected to the left and the bottom, to the top. For a gaussian window 
the matrix [ajf^] is symmetric with purely real entries. 

In equations 17 and 18 the FFTs can be used instead of 9". The 
results should be decimated according to the desired step in the 
frequency dom?.in. To achieve more accurate results a basic resolution of 
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Afa/4 has been used for the FFTs and consequently a decimation by 4:1 in 
17 and 2:1 in 18. 

The convolver block utilises the current filter coefficients and 
5 computes as its output the convolution of the reference signal, x, with 
the filter. For this purpose it is necessary to obtain the transversal 
form of the filter. The previous computation resulted in a set of 
weights. The transversal form of the filter may be obtained b/ computing 
the weighted sum of the basis functions with the recently computed 
10 weights. Tliis step may be carried out by using an IFFT followed by 

multiplication of the result with a sampled window function. Depending on 
whether the convolver works in the time domain or in the frequency domain 
it uses the computed impulse response or performs an additional FFT to 
get the filter frequency response, Q(«). 

15 

In the preferred embodiment the time domain (the lag domain) is 
used and the signals (and their correlation functions) are sa."npled with a 
period T=l/fs, Hence all expressions in the frequency domain get the form 
of a discrete Fourier transform (DFT) and the fast Fourier transform 
20 (FFT) could be efficiently used for actual computation. This yields the 
following general expression for the filter taps; 



In the preferred embodiment with a gaussian window function, the 
filter length, L, is chosen to be equal to 2N and a computation of the 
filter in the time domain is found by putting a zero between every tv;o 
coefficients Cj^ and performing a 2N-point IFFT in the equation for the 
filter taps, yielding the following expression for the filter taps ; 

= = JVc^exp{-(^)=') (c„} IFFTlCj,] (20) 

The delay on the main microphone signal, y, with to=0 is equal to N taps 
(half the filter length). 
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The equation solver block is activated every few frames and the 
system of N equations with N unknown coefficients Cj^ is solved. The 
equation solver uses an appropriate one of the many existing variations 
of Gauss elimination to solve the system efficiently in 0(N) operations. 

5 

In the case that the convolver uses the same signal frames as the 
correlation block, it may use FFTs of the reference x that have already 
been computed. The convolver can also use smaller signal frames, which, 
in order to minimise delay, is the case in the preferred embodiment. In 
10 that case the convolver should provide the computation of FFTs of smaller 
signal frames. Also the filter impulse response should be split into 
several frames and the corresponding FFT calculated. 

A filter has to be causal. Therefore in the case that t;<0 the 
15 filter is shifted by tj/fg taps. In the difference z=y-Qx, the main 
microphone signal, y should therefore be delayed as well. 

The device includes a delay 80 for the main microphone signal, 
which balances any undesirable delay in the adaptive filter due either to 
20 incomplete causality of the filter or to delay in the convolver. 

The control logic contains either short-time power or amplitude 
estimators for the reference and the subtracter output. The forgetting 
factor, a, which is utilised in the correlation block, is set to a 
25 certain fixed function of the two control variables; the short- time 

estimation of the reference, E^^\ and the short-time estimation of the 
subtracter output, E^^K 

If the main speaker signal is much larger than the noise component 
30 in the main microphone output and E^^^/E^^^ is large, a should be set to 1, 
since past segments of the signal contain information relevant for the 
current filter adaptation. If the interfering noise signal is dominant 
and E^^Ve^^^ is small, the current signal behaviour bears more relevance 
on the filter adaptation than past signal behaviour, and so a is set to a 
35 value 0<a<l. This also provides a faster adaptation. Correspondingly a 

could be chosen to be a certain function of the ratio E^^Ve^^^- However if 
both E^^^ and E^^^ are estimated to be small (i.e. both the microphones are 
at low noise level) then the set a=l is preferable, since it freezes the 
filter components, 

40 
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The adaptive filter is capable of fast adaptation because the 
adaptation time can be made as small as the filter length, 

A conventional analog AGC block, which controls the main microphone 
output before A/D conversion could be used in the device without 
additional modifications, since the adaptation rate of the proposed noise 
cancellation device is sufficient to accommodate slow gain changes, A 
digital AGC can be used in cases for which an analog AGC at the main 
microphone is absent- This digital AGC block utilises the average pov;er 
or amplitude estimate of the sub tractor output. Either a separate 
estimator could be used or the estimator in the control block could be 
utilised, possibly after additional lov;-pass (optionally non-linear) 
filtering. This non-linear filtering could provide for different time 
constants for rise-time and fall-time. The gain G is defined as a fixed 
positive function of the output E of the estimator. The larger E is , the 
smaller the gain G is, and vice versa. For small signals E=0 the gain is 
limited by a certain constant Gq. The subtractor output is thf»n 
multiplied by the gain to produce the output. 

In addition logic can be provided that digitally performs signal 
compression and decompression; i.e. the incoming signal could be 
decompressed from a non-linear (usually logarithmic) form to a linear 
form and the outgoing signal could be compressed from a linear to a non- 
linear form. 

The noise cancellation device can be used as an echo or bawling 
canceller for a full duplex speakerphone • In this application the speaker 
phone is arranged as shown in Figure 2 with the signal in from the remote 
telephone fed to loudspeaker 110 and in addition used as reference signal 
X. The signal from the microphone 100 is used as input signal y and the 
adaptive noise cancellation device 90 cancels any feedback signal from 
the loudspeaker to the microphone. 



17 
CLAIMS 

1. An adaptive noise cancellation device having a signal input (y), a 
5 . reference input (x) and a signal output (z), the device comprising: 

convolution logic (10) for convolving the signal from the reference input 
(x) with a discretised L-tap filter to form a filtered reference signal; 
and 

10 

logic (20) for subtracting the filtered reference signal from the signal 
input to form the output signal; 

logic for generating the filter taps as a linear combination of N basis 
15 functions each having a corresponding coefficient Ck; and 

logic for repeatedly determining the coefficients Ck which minimise the 
power in the output signal (z) , 

20 characterised in that 

N is less than the number of filter taps L and the basis functions have a 
portion of finite width, outside of which portion the functions are 
substantially zero, both in the frequency and time domains. 

25 

2. An adaptive noise cancellation device as claimed in claim 1 
wherein the M basis functions have the form: 

30 

for N values of an integer k, where w(t) is a window function having a 
portion of finite width, outside of which portion the function is 
substantially zero, both in the frequency and time domains. 

35 3, An adaptive noise cancellation device as claimed in claim 2 wherein 
the window function w(t) is a gaussien function. 
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4, An adaptive noise cancellation device as claimed in claim 3 wherein 
the window function has the form: 

Y 

5 with AO = 2TTfs/N and y in the range 0.25N to 0.75N- 

5. An adaptive noise cancellation device as claimed in any preceding 
claim v/herein the logic for generating the filter response function 
comprises : 

10 

correlation logic (30) for generating the cross-correlation between the 
input signal and the refernce signal and for generating the auto- 
correlation of the reference signal; 

15 logic (40) for calculating the right hand side of a system of linear 
equations for the coefficients Ck using the equation 

where G^yCa)) is the Fourier transform of product of the window function 
20 v;(t) and the cross-correlation between the input signal and the reference 
signal ; 

logic (40) for calculating the coefficients of the system of linear 
equations using the equation, 

25 

v;here G^(h,[x) is the Fourier transform of the product of the 
autocorrelation of the reference signal and an auxiliary function 

30 where W(a)) if the Fourier transform of the window function w(t); and 
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equation solving logic (50) for solving the equations ai^Cj,=rj for Cj^-. 

6. An adaptive noise cancellation device as claimed in claim 5 v/herein 
the correlation logic (30) comprises: 

5 logic for segmenting the input and reference signals into franes of a 
length equal to the filter length L; 

logic for generating a Fourier transform of the frames; 

logic for comfoining the Fourier transforms of the frames with a current 
segment of the auto-correlafcion function and cross-correlation functions 
10 v/eighted by a factor a, where a is between 0 and 1, to form an updated 
segment of the cross-correlation and auto-correlation functions. 

7. . An adaptive noise cancellation device as claimed in claim 6 
comprising control logic (60) for calculating the factor a as an 

15 increasing function of the ratio of the power or amplitude of the 
reference signal and the output signal. 

8. Adaptive noise cancellation device as claimed in any preceding 
claim wherein N is half the filter length, 

20 

9. A speakerphone capable of full-duplex operation having a signal 
input and a signal output and comprising an adaptive noise cancellation 
device (90) as claimed in any preceding claim; 

a microphone (100) arranged to provide the input signal (y) to the noise 

25 cancellation device; 

and a loudspeaker (110) connected to the signal input of the 
speakerphone, wherein the signal input of the speakerphone provides the 
reference input (x) for the noise cancellation device and the output of 
the noise cancellation device (z) provides the signal output of the 

30 speakerphone. 
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